Evaporating dynamical horizon with Hawking effect in Vaidya spacetime 
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We consider how the mass of the black hole decreases by the Hawking radiation in the Vaidya 
spacetime, using the concept of dynamical horizon equation, proposed by Ashtekar and Krishnan. 
Using the formula for the change of the dynamical horizon, we derive an equation for the mass 
incorporating the Hawking radiation. It is shown that final state is the Minkowski spacetime in our 
particular model. 
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I. INTRODUCTION 

In the study of black hole evaporation, there has 
been an important issue how black hole mass decreases 
as a back reaction of the Hawking radiation lj . We have 
to confront with this issue to resolve the information 
loss paradox. There have been many works concerning 
black hole evaporation, either in string theories |2j[3||_l(. 
or semiclassical theory typically using apparent horizon 
0| . Hiscok studied spherical model of the black hole 
evaporation using the Vaidya metric, which we also use 
in present work, to solve the black hole evaporation prob- 
lem. However, he simply set a model not taking account 
of the field equation. Hajicek's work0 treated the black 
hole mass more generally than our present case. How- 
ever, he did not use the field equation either. One of 
the more recent studies is Sorkin and Piran's work [8| on 
charged black holes. And neutral case has been done by 
Hamade and Stewart 9J. Their conclusion is that black 
hole mass decreases or increases depending on initial con- 
dition. They used a model of the double null coordinates, 
and obtained a numerical result. But they did not con- 
sider the Hawking effect directly but they used massless 
scalar field as a matter. Brevik and Haines calculated pri- 
mordial black hole evaporation^]- Very recently Hay- 
ward studied black hole evaporation and formation using 
the Vaidya metric [Tlj . It seems no analytical equation 
has been proposed for the black hole mass with the Hawk- 
ing effect taken into account. 

The dynamics for the black hole mass with the Hawking 
effect is a long standing problem. Page^^ derived the 
equation of mass intuitively, that is M oc — M~ 2 . But 
it does not come from the first principle. We will com- 
ment on his intuitive result in the final section. To derive 
the equation of mass from the first principle we should 
treat the Einstein equation with the back reaction term 
by the Hawking radiation. However, the Einstein equa- 
tion cannot be analytically solved, because the equation 
contains fourth derivative terms as back reaction. Re- 
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cently Ashtekar and Krishnan derived an equation which 
describes how the horizon changes in time. It needs only 
information of the horizon surface. 

In section ITT1 dynamical horizon is reviewed. In section 
IIHI the location of the dynamical horizon in the Vaidya 
spacetime is identified. And then in section llVl the dy- 
namical horizon equation is written down in the case of 
Vaidya matter with the Hawking effect being taken into 
account. Section [V] is devoted to conclusion and discus- 
sions. 



II. DYNAMICAL HORIZON 

Ashtekar and Krishnan considered dynamical hori- 
zon |13||14|. and derived a new equation that dictates 
how the dynamical horizon radius changes. Apparent 
horizon is a time slice of the dynamical horizon. The 
definition of dynamical horizon is, 

Definition. A smooth, three-dimensional, spacelike 
submanifold H in a space-time is said to be a dynamical 
horizon if it is foliated by preferred family of 2-spheres 
such that, on each leaf S, the expansion 0(2) of a null 
normal l a vanishes and the expansion Q( n ) of the other 
null normal n a is strictly negative. 

The requirement that one of the null expansions is 
zero comes from the intuition that black hole does not 
emit even light. And the requirement that other null 
expansion is strictly negative comes from that null 
matter goes in black holes inwards. 

In this definition, we can recapitulate the important 
formula which gives a change of the dynamical horizon 
radius by the matter flow, using 3+1 and then 2+1 
decompositions and also the Gauss-Bonct theorem. 

<§-§>= L T -^' fv 

Here Ri , R\ the radius of the dynamical horizon , T a b 
is the stress-energy tensor, \a\ 2 — a a bcr ab ', |£| 2 = CaC° 
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where a a t, is the shear, and ( a := q ab r c VJb, with two 
dimensional metric q a b, and := A^? a (Af# := \dR\) 
where R is the radius of the dynamical horizon. 
This is the dynamical horizon equation that tells how 
the horizon radius changes by the matter flow, shear 
and expansion. In the spherically symmetric case that 
we shall consider in what follows the second term of 
the right hand side vanishes. Although in the case 
of quantum field theory in curved space time, the 
dominant energy condition does not hold|lfi|||l7||. we 
can use the dynamical horizon equation because it is 
valid even when the black hole radius decreases. And 
the dynamical horizon equation is a consequence of 
the Einstein equation. We use the dynamical horizon 
equation in place of the Einstein equation. 
Because of the negative energy, we can expand the 
definition of the dynamical horizon in the case of time- 
like. Then the definition of the dynamical horizon is now, 

Definition (modified version). A smooth, three- 
dimensional, spacelike or timelike submanifold H in 
a space-time is said to be a dynamical horizon if it 
is foliated by preferred family of 2-spheres such that, 
on each leaf S, the expansion 8m of a null normal l a 
vanishes and the expansion 6( n ) of the other null normal 
n a is strictly negative. 

We can easily calculate the timelike dynamical horizon 
equation only replacing f a and f a . And the way of 3+1 
and 2+1 decomposition is replaced. The timclike dy- 
namical horizon was defined by Ashtekar and Galloway 
and called timelike membrane |15| . 



III. VAIDYA SPACETIME 

The Vaidya metric is of the form 

ds 2 = -Fdv 2 + 2Gdvdr + r 2 dfl 2 , (2) 

where F and G are functions of v and r, and v a is null 
vector and r is the area radius, and M is the mass defined 
by M = §(1 — (fs)) a function of v and r. This metric is 
spherically symmetric. By substituting the Vaidya met- 
ric (2) into the Einstein equation so that we can identify 
the energy-momentum tensor T a b as 



8ttT vv := -^(FM, r + GM tV ) 



&TtT rr 
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(3) 
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We do not need to check that the solution of the dy- 
namical horizon equation satisfies the Einstein equation. 
Because we would like to consider the Schwarzschild like 
metric, we set v = t + r* , where r* is tortoise coordinate 




FIG. 1: For the case that the dynamical horizon decreases, 
we should choose l a = — t a + f a so that l a points into the 
dynamical horizon. 



with dynamics 



r* = r + 2M (v) In 



2M(v) 

For later convenience, we write, 

dr 



- 1 



dr* 



(6) 
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There are two null vectors, 



l a = 



r 
i' 



(8) 



corresponding to the null vector v a , and the other is 





( 


n r * 




n e 









F a- 1 



F — 2Ga 






(9) 



/ 



Here we multiply a 1 so that l a = v a . This choice of the 
null vector l a is explained in figure ^ From now on we 
put, 



G = 1, 



2M{v) 



(10) 
(11) 



in a similar form to the Schwarzschild metric, assuming 
that M(v) is a function of v only. For a constant M, 
the metric coincides with the Schwarzschild metric. We 
calculate the expansions 8ni and 6( n ) of the two null 
vectors l a ,n a , because the definition of the dynamical 
horizon requires one of the null expansions to be zero 
and the other to be minus. The result is, 



e 



(0 



B( n) = 



-(2F-a) 
r 

1 f-2F 2 + aF-2a 2 
~a~r V -F + 2a 



From 0(/) =0 we get, 



2F- a = 0. 



(12) 
(13) 

(14) 
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we can check that the other null expansion 0(„)is strictly 
negative. Therefore in this case, we can apply the dynam- 
ical horizon equation. In the usual Schwarzschild metric 
with dynamics, both expansions become zero. This is the 
one of the reasons why we choose the Vaidya metric. By 
inserting equation (7) to equation (14), we obtain 



a = r| 1 - 2M„ln - 1 
' \2M 



-M, 



M(r/2M - 1 



(15) 



Note that a is proportional to F. Now we solve = 0, 
to determine the dynamical horizon radius as 

2F-a = 2F 

Fi 1 - 2M„ln (— 1 

' v \2M 



M{r/2M - 1) 
= 0. 



(16) 



From this equation we obtain, 



I + I - 2M v In 

I'D 



M(r D /2M - 1) 



\2M J 



M v = 0. 



(17) 



Here F = is also the solution of the dynamical horizon. 
The dynamical horizon radius r\o is given by solving (17) 
as 



r D = 2M + 2Me~ v/2M . 



(18) 



Note that this dynamical horizon radius is outside the 
r = 2M. that is othor solution. 



IV. DYNAMICAL HORIZON EQUATION WITH 
HAWKING RADIATION 



At first, we should derive the energy- momentum 
tensor T|, for the integration of the dynamical horizon 
equation. For this end we derive it from the given Vaidya 
matter. For G — 1, F = 1 — 2M ( V > : the non- vanishing 
components of the energy-momentum tensor becomes 



T — 

Tlr* = 

T , « = 



1 



(FM, r + M tV ) 



Anr 2 
0. 



(19) 

(20) 
(21) 



Here we have made the coordinate transformation from 
r to r* . Writing T t i in terms of T vv and T vr * given by 



(PUm with 1° 



T t i = -T, 
1 



v , we see 

(-FM r - M„ - aM tr ) 



4 7rr 2 

47rr 2 2 



(22) 



With i a being the unit vector in the direction of t a , we 
obtain 



Til 



1 5 



-M.„F- 



Anr 2 2 

For the dynamical horizon integration l|T|l. we get 



(23) 



^ r D T tl dr D 

where we have used 

F = 

and the fact 



M 2 



' (l + e- v / 2M )dM, (24) 



Mi 



-v/2M 



dM 
dv 



-v/2M 



2(1 



-v/2M ■ 



-v/2M\ 



V 



J_ e -v/2M 



(25) 



,(26) 



changing the integration variable from to M. In the 
above calculation, we treat M tV and F -1 with ro fixed, 
because these functions are used only in the integration. 
Inserting equation (|24() to the dynamical horizon equa- 
tion JJJ, we obtain 



-(2M + 2Me~ v/2M ) 



M 2 

Mi 
M 2 K 

Ul + e-"/ 2M )dM. 
Mi 2 



Taking the limit M2 
3 



Mi = M, we obtain 



2 (l + e -^M ) + _i_ e -/ 2 M =0 _ 



(27) 



(28) 



This equation is the dynamical horizon equation in the 
case that only the Vaidya matter is present. There is no 
solution of this equation, except the trivial one (F = 
or r = 2M), so 



r D = 2M(v). 



(29) 



Here M(v) is the arbitrary function only of the v, which 
represent the Vaidya black hole spacetime. 
Next, we take into account the Hawking radiation. To 
solve this problem, we use two ideas that is to use the 
dynamical horizon equation, and to use the Vaidya met- 
ric. The reason to use the dynamical horizon equation 
comes from the fact that we need only information of 
matter near horizon, without solving the full Einstein 
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equation with back reaction being the fourth order dif- 
ferential equations, for a massless scalar field. For the 
matter on the dynamical horizon, we use the result of 
Candelas [l^, which assumes that spacetime is almost 
static and is valid near the horizon, r ~ 2M. 



T-a — — Tu 



1 



dtouj 3 



2tt 2 (1 - 2M/r) J e 87rM "> - 1 
1 



2cM 4 tt 2 (1 - 2M/rY 
where we have used a well known result, 



(30) 



1 



15a 4 



(31) 



and where c = 61440. This matter energy is negative 
near the event horizon. In the dynamical horizon equa- 
tion, if black hole absorbs negative energy, black hole 
radius decreases. This is one of the motivations to use 
the negative energy tensor. Next we replace length of t to 
unit length, because in the dynamical horizon equation t 
is used, so 



and therefore, the energy tensor becomes 
T- 1 



2M 4 ctt 2 (1 - 2M/r) 2 



(32) 



(33) 



Calculating the integration on the right hand side of Q 
for this matter, 



' D 



= b 



M 4 (l - 2M/r D ) 2 
4Af 2 (l + e- l '/ 2M ) 4 dr D 



M i e -v/M 



dM 



dr D 
dM 



4(1 + e -„/2Af )4 ^ /M 



M 2 



X (2(1 + e~ v ' 2M ) + ^e-^ 2M )dM. 



(34) 



Here we insert the expression for rjj (|18(l in the first 
line, and the expression for dro/dM = 2(1 + e~ v / 2M ) 
v_ e -v/2M j g usec ^ Here b is a constant calculated in [l 



6 = 



1 



30720tt 



(35) 



If we also take account of the contribution of the Vaidya 
matter, and inserting this into the integration to the dy- 



namical horizon equation we obtain 
-(2M + 2Me~ v/2M ) 



M 2 



M 1 



M 2 n2 



2 2 (l + e- /2M )V ,/M 



Mi 



M 2 



x ^2(1 + e - v ' 2M ) + ^ v,2M ^j dM 

+ Ul + e- v / 2M )dM. 
J Mi 2 

Taking the limit M 2 -» M x = M, we finally get 



(36) 



[ {l+e -v/2M ) + 



u -v/2M 
2M 

2 2 (1 + e-"/ 2M ) 4 
M 2 



xe"/ y ^2(1 + e - v ' %M ) + Yi e ~ V ' 2M ^j ' (37) 



or 



M 2 = 



{l + e -v/2M^ e v/M 



i( 1 + e -v/2M} + v e - v/ 2M 



X (1 + e 



-v/2M\ 



,-v/2M 



2M 



(38) 



This is the main result of the present work that describes 
how the mass of black hole decreases. This equation 
is the transcendental equation, so usually it cannot be 
solved analytically. However, with the right hand side 
depending only on —v/2M, we can easily treat Ea. l|38|) 
numerically. Figure [21 is a graph of M as a function of v 
If the dynamical horizon were inside the event horizon, 
the dynamical horizon radius would be 

r D = 2M - 2Me~ v/2M . 

In this case, the dynamical horizon equation would be- 
come 



M 2 = 



-8b(l-e- v / 2M ) 4 e v / M 
|(1 - e-/2Af) _ _2_ e -«/2M 



(I _ e -v/2M\ v _ -v/2M 

V ' 2M 



The singular behavior of this expression excludes its 
physical relevance. 

Now we show an approximation of the Ea. (|38|l in par- 
ticular limiting case. Taking the limit M 0, and 
—v/2M = const, we can see that (|5*51l becomes, 



bCi 
M 2 



bvC 2 



= 0. 



Where C\, C 2 are positive constants, or 



(39) 



(40) 
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FIG. 2: Numerical calculation of the black hole mass M as a 
function of v from the equation (60) 



So, in the vanishing process the mass is proportional to 
v. For M — > large 



M 



C, 



logM' 



(41) 



where C3 is a positive constant. It comes from the limit 
M —> 00 and —v/2M — > 00. In this limit the equation 
(|38|l becomes 1; = — 2MlogM. This is different from 
Page's result [l2|- Because if M goes to large, the dy- 
namical horizon radius increases as M 2 , so absorbed en- 
ergy also become large. From this reason derivative of 
M by v changes. If we do not consider next order, the 
derivative of M becomes M — — C4, so that Airr^T 4 « 1, 
contradicting with Page's intuition. 



From the analysis of the transcendental equation (38), we 
have shown that the black hole mass eventually vanishes 
and the spacetime becomes the Minkowski spacetime in- 
dependent of the initial black hole mass size. 
The dynamical horizon method in this paper can take 
into account of the back reaction of the Hawking radi- 
ation without solving the field equation which contains 
the fourth order differentials. 

In the limit of the black hole mass going to zero, the 
derivative of the mass becomes small in proportion to 
the null coordinate (v = t + r*). On the other hand 
as the black hole mass becomes large, the derivative be- 
haves the minus of the inverse of the logarithm of the 
mass. Our result, which is different from Page's result, 
comes from the fact that in the large mass limit, the black 
hole radius behaves like quadratic of the black hole mass. 
This probably comes from when large mass limit that the 
approximation r — > 2M is broken. 

We would like to compare the present work to the pre- 
ceding works. Sorkin and Piran or Hamade and Stewart 
used a massless scalar field instead of the Hawking ra- 
diation as the back reaction directly. The conclusion of 
their paper is that black hole starts with the small mass 
and it evaporates or increases. However, it is shown in 
the present work that even if the black hole starts with 
a large mass it always vanishes. 

Although we have treated the black hole evaporation 
semi-classically, we hope this work will give an intuition 
to quantization of black holes. 



CONCLUSION AND DISCUSSIONS 



We have derived an equation which describes how 
the black hole mass changes taking into account of the 
Hawking radiation, in the special Vaidya spacetime which 
becomes the Schwarzschild spacetime in the static case. 
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